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Abstract 

A diffraction theory in a system consisting of left-handed and right-handed materials 



> 

£T) ' is proposed. The theory is based upon the Huygens's principle and the Kirchhoff 's 

integral and it is valid if the wavelength is smaller than any relevant length of the 
system. The theory is applied to the calculation of the smearing of the foci of the 
fS| ' Veselago lens due to the finite wavelength. We show that the Veselago lens is a 

unique optical instrument for the 3D imaging, but it is not a "superlens" as it has 
, been claimed recently. 

a. 
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In his seminal work Veselago [1,2] has introduced the concept of left-handed 
^ ■ materials (LHM's). In a simplest case the LHM's are materials with simulta- 

neously negative electric permittivity e and magnetic permeability fi in some 
frequency range. In the LHM the vectors k, E, H form a left-handed set, while 
in the usual materials (e > 0, fi > 0) they form a right-handed set. If imagi- 
nary parts of e and fi are small, the electromagnetic waves (EMW's) propagate 
in the LHM but they have some unusual properties. All these properties orig- 
inate from the fact that in the isotropic LHM the Poynting vector S = E x H 
is anti-parallel to the wave vector k. 

Consider a propagation of the EMW from a point source located at the point 
z = —a through an infinite slab of the LHM with the thickness d and a usual 
right-handed material (RHM) at z < and z > d (Fig.l). It is obvious that 
S z > everywhere at z > —a because the energy propagates from its source. 
The directions of vector k for different rays are shown by arrows. They should 
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be chosen in such a way that at both interfaces tangential components of 
vector k for incident, reflected and refracted waves are the same. Another 
condition is that the component k z should be parallel to S z in the RHM 
and anti-parallel in the LHM. Then in the LHM k z is negative. It follows 
that the Snell's law for the RHM-LHM interfaces has an anomalous form: 
smi/ sinr = —n'/n, where i and r are the angles of incidence and refraction 

respectively, n! = ^\e'\\fi'\/eoHo and n = yjefx/ eo/^o are positive refractive 
indexes for LHM and RHM respectively. The angles of reflection are equal to 
the angles of incidence. 

The device shown at Fig. 1(b) is a unique optical lens proposed by Veselago. 
In this lens e = — e' and fi = — //, then n' — n and % = —r. It is easy to show 
that in this case the reflected wave is completely absent. Since all the rays 
going right from the source have % = — r, all of them have foci at points z = a 
and z = 2d — a as shown in Fig. 1(b). 

All the ideas above have been put forward by Veselago about 35 years ago [1,2]. 
Recently the method of fabricating of the LHM's on the basis of metallic 
photonic crystals has been found and the San Diego group has reported the 
first observation of the anomalous transmission [3] and even the anomalous 
Snell's law [4]. Both observations have been interpreted as a result of negative 
e and ji. 

One can see that the Veselago lens, shown at Fig. 1(b), is an absolute instru- 
ment because it images stigmatically a three-dimensional domain — d < z < 
and the optical length of any curve in the object space is equal to the optical 
length of its image [5]. Note, that the definition of the absolute instrument 
assumes geometrical optics only. Since the LHM's have been already obtained 
we think that the Veselago lens can be extremely important device for the 3D 
imaging. 

Pendry [6] claims that the Veselago lens has another unique property. Due to 
Pendry the resolution of this lens does not have a traditional diffraction limi- 
tation which follows from the uncertainty principle. Pendry has introduced a 
new term "superlenses" with the Veselago lens as a first representative of this 
class. Several works have recently appeared where the concept of superlenses 
was criticized [7-10]. On the other hand, Ziolkowski and Heyman [11] proposed 
a refined version of Pendry 's analysis which confirms his result. Haldane [12] 
interpreted results of Pendry in terms of surface waves described by Ruppin 
[13]. Haldane also proposed a way how to overcome some mathematical incon- 
sistencies of papers [6,11], shown before in [14,10]. Haldane claims that even 
though the perfect focusing is impossible the diffraction limit can be defeated. 

In this paper we propose a general scalar theory of diffraction in the LHM 
which is based upon the Huygens's principle and the Kirchhoff's integral. As 
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any diffraction approach our theory works under condition that the wavelength 
is much smaller than any relevant geometrical length in the problem. We 
apply this theory to the Veselago lens and calculate the smearing of the foci 
due to the finite wavelength. Thus, our result does not support the idea of 
"superlens" . 

The Green function for the Helmholtz equation for the LHM, which describes 
propagation of a spherical wave from a point source, has a form exp {—ikR)/R, 
where k = um/c and R is a distance from the source. At a small element of the 
sphere R = const the spherical wave can be considered as a plane wave which 
is characterized by the Poynting vector S and wave vector k both with the 
radial component only. Since S is directed along the external normal to the 
surface element, the wave vector k in the LHM is directed along the internal 
normal. It is easy to see that our Green function obeys these properties. 

Following the principles of the scalar theory of diffraction [15] the field u at 
the observation point P can be written in a form of a surface integral 



where R is the length of the vector from the point P to the surface element, 
df n is the projection of the surface element df on the plane perpendicular to 
the direction of the ray coming from the source to df, b\ is a constant for any 
LHM. To find b\ one can consider a plane wave with the wave vector normal to 
the infinite plane of integration. Since this plane is fictional, the constant can 
be found from the condition that the Huygens's principle in the form Eq.(l) 
reproduces the same plane wave. Doing the calculations similar to Ref.[15] one 
gets bi = —k/2ni, so that for the LHM the constant bi has a different sign 
than the similar constant b r for the RHM. 

The Huygens's principle can be applied to any interface which has a curvature 
larger than the wavelength. It gives the correct direction of refracted waves 
but it does not give the amplitudes of both refracted and reflected waves. 
However, it can be successfully applied to the Veselago lens where reflected 
waves are absent. 

Note, that there are other methods to describe the diffraction which may 
be also used if the source of the rays is unknown. They are described and 
compared in details in Jackson's textbook [16]. One can show that all the 
methods give the same result at r = —i. 

Moreover, the Kirchhoff integral Eq. (1) may be considered as a hint to get a 
solution, valid at large k. The function u as given by Eq. (1) is meaningful only 
at the distances from the interface larger than the wavelength. Nevertheless, 
if we consider solution obtained for these distances exactly at the interface, 




ikR 



(1) 
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we find that it satisfies proper boundary conditions. Namely, u is continuous 
and du/dz changes sign at the interface (See[ll]). 



Now we apply Eq.(l) to the Veselago lens. To find the field u inside the slab 
we shall integrate in Eq.(l) over the plane z — 0. The field u in this plane is 
produced by a point source and has a form 



u(x,y,0) 



gik^J a 2 +x 2 +y 2 

\J d 1 + x 2 + y 2 



(2) 



The field inside the slab can be found using Eq.(l) with a constant b r i instead 
of bi because now we are integrating over the RHM-LHM interface rather 
than over the fictional surface in the LHM. In a similar way at the LHM- 
RHM interface one should use a constant bi r . Using the method described in 
Ref.[15] it is easy to show that b rt = b t and b lr = b r . Thus one gets 



u(x,y,z) = ha J J ^ 



iky/ 'a 2 +x 2 +y 2 -ik^J z 2 + (x 1 -x) 2 +(y±-y) 2 

+ + Vi ^/z 2 + (x 1 -x) 2 + (y 1 - y y 



dxidyx, (3) 



where the additional factor a/ \J a 2 + x\ + y\ is the cosine of the angle between 
the ray, coming from the source to the point {xi,yi,0} and the unit vector 
in z direction. One can see that the optical lengths for all rays (the sum of 
exponents in the integrand of Eq.(3)) from the point source to the focus, 
located at z = a, x = y = 0, are zero and the value of the field at the focal 
point u(0, 0, a) = ik, while the geometrical optics gives an infinite field in this 
point. To find u in the vicinity of the focus one should expand the integrand 
in Eq.(3) near the point (0,0, a) assuming x a, y a and |£| a, where 
C = z — a. Then one gets 



ka 



2-7T OO 



2ni 



exp 



ik 



pr cos (ip) — a( 



Vd^T 







( a 2 + r 2)3/2 



rdipdr, 



(4) 



where p 2 = x 2 + y 2 . Equation (4) gives 



u(p,C) = k 



i sin {k\J p 2 + ( 2 



— J Jo(kpVi — s 2 ) sm(k(s)ds 



(5) 



At p = one gets an analytical expression 

1 — cos (k() + i sin (k() 



«(o,C) = 



c 



(6) 



Another analytical expression can be obtained at ( = 0: u(p, 0) = isin(kp) / p. 
The same result has been obtained in Ref. [10]. Figure 2 shows dimensionless 
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function \u(p, ()\ 2 /k 2 as given by Eq.(5). One can see that the smearing of 
the focus is anisotropic. The half-width in z direction is approximately one 
wavelength while in p direction it is approximately twice as less. 

Now we find the field u in the close vicinity of the second focus located at 
x — y — 0,z — 2d — a. The general expression for u at z > d differs from 
Eq.(3). One should apply the Huygens's principle to both interfaces located 
at z = and z = d. Then the expression for the field has a form 

u(x,y,z) =bib r ad J J dx x dy x J J dx 2 dy 2 - 



-oo — oo 



a z + x\ + yf 



-iky/d 2 +(x 1 -x 2 ) 2 +(yi-y2) 2 -ife^/ (z-d) 2 +(x2-x) 2 + (y2-y) 2 



d 2 + ( Xl - x 2 ) 2 + (yi - y 2 ) 2 ^_ d )2 + ( a . 2 _ a .) 2 + ( y2 _ y )i 



(7) 



To calculate these integrals in the vicinity of the second focus using inequalities 
kd ^> 1, ka >> 1 one should introduce new variables {s, t} instead of variables 
{ X 2,U2} by relations 

x 2 = -(--l)x 1 + s, y 2 = - l) yi + t. (8) 
a a 

One can see from Eq.(7) that at s = t = optical lengths of all rays exiting 
from the point source at z — —a and coming to the second focus at z = 2d — a 
are equal to zero. Thus the new variables {s,t} describe deviation from the 
geometrical optics and they should be small. Introducing the new variables 
and expanding the exponents in Eq.(7) one can get an expression for the field 
u in the vicinity of the second focus. For r\ = z — 2d + a, \rj\ a one gets 

u(x, y, 7]) = -u(x, y, C)*le=r?> ( 9 ) 

where the function u(x,y,Q is given by Eq.(3). In other words, the smearing 
of the \u(x, y, i])\ 2 in the second focus is the same as the smearing in the first 
one. 



Our result has a simple interpretation in terms of uncertainty principle. One 
can see that the main contributions to the integrals in Eqs.(3,7) come from 
a circle of a radius a in the x-y plane around z axis. Thus, the effective 
aperture of the Veselago lens is of the order of a. The uncertainty of the 
lateral momentum is of the order of 1/a, which gives the smearing of the foci 
of the order of 1/k . 

Now we compare our results with the analytical calculations of Pendry [6] 
and Ziolkowski and Heyman [11]. Our paper is valid for large k only while in 
[6,11] the near field region is also considered. The main difference between our 
works is that in Ref. [11] the incident spherical wave is expanded into series of 
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Fig. 1. Reflection and refraction of light outgoing from a point source at z = —a and 
passing through the slab of the LHM at < z < d. Refraction of light is described 
by the anomalous Snell's law. The arrows represent the direction of the wave vector. 
The reflected waves are shown by dashed lines near each interface only. The slab 
is surrounded by the usual RHM. (a) n' > n. (b) The Veselago lens. The reflected 
waves are absent, all rays pass through two foci. 

the plane waves in the x — y plane. This series contains the evanescent waves 
(EW's) which are claimed to provide the "superlensing" . The matching of 
these plane waves at the interfaces is done independently. Note that each EW 
has the Poynting vector in the plane of the lens, which does not correspond to 
the problem under study. Moreover, the superposition of the EW's diverges 
within three-dimensional domain a < z < 2d — a [10,12,14] at large values 
of k. This happen because their amplitudes increase as exp(A;,2) with positive 
k. Haldane[12] argues that the "ultra violet cutoff" for this problem exists, 
but it is outside macroscopic electrodynamics, so that upper value of k is very 
large. However, even with the cutoff the solution looks strange. The field is 
exponentially large within the domain a < z < 2d — a which occupies all 
the space between the two foci. The field in each focus has a different limit 
from the left and from the right side. Therefore, we do not think that such a 
solution means a perfect focusing. 

Finally, we have proposed the theory of diffraction in a system, consisting of 
the LHM and the RHM and have applied this theory to the calculation of the 
smearing of the foci of the Veselago lens. This smearing is of the order of the 
wavelength. 
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